Resonant tunneling through a macroscopic charge state in a 
superconducting SET transistor 

D.V. Averin 

Department of Physics, SUNY at Stony Brook, Stony Brook, NY 11794, USA 

A.N. Korotkov 

Institute of Nuclear Physics, Moscow State University, 119899 GSP Moscow, Russia 

A.J. Manninen and J. P. Pekola 
Department of Physics, University of Jyvaskyla, P.O. Box 35, 40351 Jyvaskyla, Finland 

Abstract 

We predict theoretically and observe in experiment that the differential con- 
ductance of a superconducting SET transistor exhibits a peak which is a 
complete analogue in a macroscopic system of a standard resonant tunneling 
peak associated with tunneling through a single quantum state. In particular, 
in a symmetric transistor, the peak height is universal and equal to e 2 /2irh. 
Away from the resonance we clearly observe the co-tunneling current which in 
contrast to the normal-metal transistor varies linearly with the bias voltage. 
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Charging effects in systems of small Josephson junctions are quite well understood by 
now - see, e.g., []1]|§. Interest, however, has been focused mostly on the interplay between 
the charging effects and Cooper pair transport, which can be described in generic terms as 
the quantum dynamics of the Josephson phase difference. The aim of this work is to study 
the quasiparticle transport in a superconducting SET transistor - a system of two junctions 
connected in series (see inset in Fig. 1). We show that the BCS singularity in the density 
of states of superconducting electrodes of the junctions brings about several interesting new 
features of quasiparticle transport. Most notably, in the vicinity of the threshold voltage Vt 
for classical tunneling the quasiparticle transport is identical to resonant tunneling through 
a single macroscopic quantum state of the transistor. 

In this work, we study the low voltage regime, V < V t , where the quasiparticles do not 
have enough energy to enter the central electrode of the transistor and can traverse it only by 
quantum tunneling through the energy barrier created by the charging energy of the central 
electrode. The effects of the superconducting density of states in the classical sequential 
tunneling were discussed recently in Ref . H . The dominant contribution to the current I in 
the regime of quantum tunneling comes from the so-called inelastic co-tunneling, the process 
in which two different electrons tunnel simultaneously in the two junctions of the transistor, 
and can be written f|] as: 

/ = e(r(V) - V(-V)) , 

IW = ^ / nj =1 [dejN^l - fie,))} 5(eV - £ e,) 
a j j=1 
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Here ex, e 2 and e 4 are the energies of the states between which electrons are transferred in 

the first and second junction, respectively, T 12 and Nj(ej) are the corresponding tunneling 

amplitudes and densities of states in the electrodes, /(e) is the Fermi distribution function, 

V is the bias voltage, and we have assumed for simplicity that all electrodes have the same 

energy gap A. 
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We restrict our attention to the case of low temperatures, T A, when the non- 
vanishing quasiparticle current exists only at large voltages, V > 4A/e, sufficient for the 
creation of quasiparticles in the two junctions. In this voltage range the energies E\£ of the 
intermediate charge states in eq. ([!]) are: 

Bk=Ec- KeV - 4A) - ^ , E 2 = E c - (1 - A) (eV - 4A) + ^ , (2) 

t^E 

where E c = e 2 /2Cs with = C\ + C 2 + C 9 denoting the total capacitance of the central 
electrode of the transistor, A = (C2 + C g ) /Cs gives the fraction of the bias voltage that drops 
across the first junction, and Qo = e x {V g C g /e + A(2A — 1)/ Ec} with {x} = x — [x + 1/2] 
can be interpreted as the charge induced by the gate voltage V g into the central electrode. 

Integrating over ei, e 2 at fixed e% + e 2 and similarly over €3, €4 we can express the co- 
tunneling rate (P in terms of the "seed" I — V characteristics Ij(U), j = 1,2, of the two 
junctions at a fixed voltage U across a single junction and no charging effects: 

r(V) = [ de ^ UV-e/e) 

K ' 27ie 2 J 1 - exp(-e/T) 1 - exp(-(eV - e)/T) 1 1 ' 

M = ^ + e - 2A + £ 2 + - e - 2A ' ^ 
^From eq. (|) we see directly that a jump of the quasiparticle current Ij(U) at U = 2A/e 
in superconducting junctions changes the voltage dependence of the co-tunneling current for 
V close to 4A/e from cubic (r(V) oc V 3 for a normal-metal transistor Q) to linear. Indeed, 
for T<A we can approximate I(U) near the threshold U = 2A/e as (see, e.g., ||): 

I(U)= I 5 Q{U - 2A/e) , /; = ^|-, (4) 

where i?j is the normal-state tunnel resistance of the jth junction. Equations (|3]) and 
give for low temperatures and eV — 4 A <C A, E c : 

, { V) = enV) = ^l(± + ±nV-^)- (5) 

Z7T hi\ h/2 e 

When the bias voltage approaches the threshold V* of classical sequential tunneling, where 
one of the energy barriers Ej vanishes, the co-tunneling current grows and crosses over into 



the current carried by sequential tunneling, in which quasiparticles traverse the transistor 
by two independent jumps across the two junctions. It is known that the energy width of 
the crossover region between the co-tunneling and sequential tunneling is determined by 
the lifetime broadening of the intermediate charge states E\$ |6|-|9|]. If the gate voltage 
is not close to the special point Qq = (1/2 — A)e where E\^ vanish simultaneously (the 
situation that corresponds to the maximum threshold voltage Vt = (AA + 2Ec)/e), then the 
current through one intermediate state, for instance E\ = E, dominates near the tunneling 
threshold. The current in the transition region can be described in this situation by simply 
adding the lifetime broadening 7 of the intermediate state in eq. (|3]) for the co-tunneling 
rate 



1 ^r T / e x,/Cx T ,eV — e. _ ,eV — e s , ,„ s 

M = - t , 7 = — A - coth — +/ 2 coth — — . 6 

E + e-2A + i7' ; 2e L V K 2T J y e ' v 2T n KJ 

(This simple approach neglects only the renormalization of E and 7 significant at tempera- 
tures exponentially small on the scale of Eq 

Combining eqs. (||), @, (f|), and (|6]) we can calculate the differential conductance of the 
transistor at low temperatures: 



G dl hhh 



dV 2tt 



A 1- A 

+ 



(7) 



_(E - X(eV - 4A)) 2 + 5 2 (E + (1 - X)(eV - 4A)) 2 + 5 2 
where E = (e/2 — \Qo\)e/Cs is the Coulomb energy barrier at V — 4A/e, and 5 = fi{I\ + 
I-i)j1e is the energy width of the charge state due to tunneling. If we use the second 
equation in the expression (f|) we see that 5 = nh^R^ 1 + i?^" 1 )A/4e 2 . Since the ideology of 
co-tunneling is applicable only to junctions with small tunnel conductance, R" 1 e 2 /h, this 
means that the width of the charge state is small, (5< A, and eq. (^) describes the narrow 
conductance peak located at the threshold V t of classical tunneling (eV t = 4A + E / X). This 
peak corresponds to the rapid current rise from almost zero to /1-Z2/ {h +I2) at V — Vt. The 
maximum conductance is achieved when E = (i.e., when the tunneling threshold reaches 
minimum) and V = Vt = 4A/e: 

dV 2nh (h + I 2 ^ 2 
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Equation (§) shows that in a symmetric transistor, where Ji = I2, the differential conduc- 
tance reaches the absolute maximum e 2 /2ith which is independent of A, Ec, or the junction 
resistance R. This universality is similar to that of the resonant tunneling through a single 
microscopic quantum state, and is quite remarkable in view of the fact that in the present 
context the quantum state is the macroscopic charge state of the central electrode of the 
transistor. 

If the energy barrier Eq is large on the scale of the width 5 of the charge state, 5 starts 
to increase with increasing Eq, i.e. increasing threshold voltage Vt- The conductance peak 
can be described analytically in this regime by retaining only the first, resonant, term in eq. 
(0), and taking into account that the peak width 5 depends then on its position V t through 
the dependence on V t of the contribution of the current J 2 through the second junction to 
S: h = h(V t -2A/e). 

The shape of the conductance peak in a symmetric transistor (with R± = R2, and 
A = 1/2) calculated numerically from the eqs. (|2|), (|j), and (|6|) without the approximation 
(||) or restrictions on E Q is shown in Fig. 1. We see that this, more accurate, calculation 
preserves all the qualitative features of the simple analytical expression (|7|): maximum 
conductance is e 2 /2nh when E = and decreases to approximately half this value at 
nonzero Eq. 

For the results discussed above to be valid, the lifetime broadening of the resonant 
charge state should not only be much smaller than the superconducting gap A, but also 
much smaller than the typical energy distance (on the order of Eq) to the excited charge 
states of the central electrode of the transistor. The condition for this is: 

A 7ih , , , . 

a = — — (i?r 1 + i? 2 - 1 )«l. (9) 

h, c c 

It is important that this condition can be violated not only when the normal-state junction 
conductances are large, but also when the energy gap A is large in comparison to the 
charging energy Ec- If it is indeed violated, the charging effects are eventually washed 
out by quantum fluctuations and the current rise at V = 4A/e becomes infinitely sharp 
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(provided that the singularity of the density of states at the energy gap A is not smeared 
out by some internal mechanism like spin- flip scattering). 

To test these predictions experimentally we fabricated and measured four superconduct- 
ing SET transistors with differing parameters. The transistors were fabricated by electron 
beam lithography on oxidised silicon by the standard shadow evaporation technique using 
aluminum electrodes and aluminum oxide junction barriers. The geometry of each of the 
transistors was such that the length of the central island was 1 jttm, its width was 80 - 120 
nm, and the overlap at the two ends of the island with the external electrodes was nominally 
70 nm. The gate electrode was simply a 100 nm wide finger, pointing orthogonally to the 
center of the island at a distance of about 0.5 /im.The gate capacitance was about 0.02 fF. 

Tunnel resistance R of the transistor junctions was measured from the large-voltage 
asymptote of the I — V characteristic of the transistor assuming equal resistances of the 
two junctions. Although we did not carry out any systematic study of how symmetric 
the transistors were, we checked from the gate voltage dependence of the threshold voltage 
V t that sample 1 had equal parameters to within 30%, and we do not expect the other 
transistors to be worse in this respect since their dimensions were larger than in sample 1. 
The charging energy Ec was measured as a half of the amplitude of the V t modulation by 
the gate voltage, and A can be obtained from the onset of the current at 4A/e. All these 
parameters of the four samples are shown in Table 1, together with the combined parameter 
a defined in eq. (^j as a small parameter of the present theory. 

Measured I — V characteristics and traces of the differential conductance of sample 1 as 
a function of the bias voltage V are shown in Fig. 2 for several values of the gate voltage. 
The curves agree qualitatively with the predictions of the theory described above. The 
differential conductance has a narrow peak of the roughly correct width at the threshold 
of classical tunneling. The height of the peak away from the resonance is slightly below 
one-half of e 2 /2nh. The main discrepancy between the experimental results (Fig. 2b) and 
the simple model calculations (Fig. 1) is that at resonance the conductance does not reach 
the ideal maximum value e 2 /2irh but rather is about one half of this value. Although 
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the asymmetry of junction resistances contributes according to eq. (j8|) to suppression of the 
resonance conductance, the actual asymmetry of our transistors was too small to account for 
the observed magnitude of this suppression. This discrepancy can be qualitatively explained 
by the fact that in close similarity to regular resonant tunneling, the resonant tunneling 
through the macroscopic charge state discussed in this work is very sensitive to all sources of 
inelastic scattering. For instance, we checked numerically that weak spin-flip scattering with 
the rate r^ 1 = O.OlA/ft is sufficient to suppress the resonant conductance peak to the level 
found in the experiment. The model with spin-flip scattering, however, did not reproduce 
correctly the full shape of the observed conductance curves and we think at the moment that 
in our transistor the resonance is suppressed by a combination of several inelastic scattering 
mechanisms including fluctuations of the bias and gate voltage (associated with the finite 
impedance of the voltage leads), and inelastic tunneling through the tunnel barriers. We 
could not characterise quantitatively all these sources of inelastic scattering, and therefore 
did not attempt to find a theoretical fit to the curves in Fig. 2. 

The results of measurements for all four samples are summarized in the right-hand-part 
of Table 1, which shows two characteristic values of the differential conductance in units of 
e 2 /2nh\ (i) Gq )£XP1 the conductance at bias voltage just above 4A/e and at gate voltage that 
corresponds to the maximum threshold voltage V t , and (ii) Gi )exp , the peak conductance 
at resonance (when V t reaches minimum). Variation of the peak conductance G\ with the 
tunnel resistance R and charging energy Ec described by Table 1 confirms that when the 
relative width of the charge states of the transistor (characterised by the parameter a of 
eq. (§)) becomes considerable, G\ increases gradually beyond e 2 /2nh. At large a, when 
G\ is much larger than e 2 /27rh, the charging effects are completely washed out by the 
quantum fluctuations of charge on the central electrode of the transistor and the differential 
conductance becomes insensitive to the gate voltage. This case is approached by sample 4 
with the largest a in which V t is practically independent of the gate voltage, and we could 
assign only one value of the characteristic conductance to this sample. 

When a is small and the charging effects are well-pronounced, the threshold conductance 
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Go originates only from the process of co-tunneling, and is much smaller than the peak con- 
ductance G\. It can be calculated from eq. (|5[) which predicts that at V = 4A/e the 
co-tunneling conductance of the superconducting SET transistor should increase abruptly 
to a finite, voltage-independent level which also does not depend on temperature at low 
temperatures. This behavior is indeed found in our three samples with larger tunnel re- 
sistances. Figure 3 shows for example the data for sample 2. At gate voltages which 
correspond to the thresholds Vt close to maximum we see the kink in the I — V curves and 
the step in the dl/dV at the onset of the quasiparticle current at V = 4A/e. (For other 
values of the gate voltage small current peaks due to Cooper pair tunneling that are visible 
in Figs. 2a and 3, overlap with the onset of quasiparticle current and do not allow to iden- 
tify the conductance jump.) The data shown in Figs. 2 and 3 were taken at temperature 
of about 100 niK. We checked that the jump in the quasiparticle conductance is practically 
temperature-independent for temperatures up to 0.4 K. 

Table 1 contains a comparison between the observed co-tunneling conductance G 0iSXp and 
Go,theory calculated from eq. (§) under the assumption of a symmetric transistor. Taking 
into account that any asymmetry of the junction tunnel resistances increases Go we can say 
that the agreement between Gq^ xp and G jheary is reasonable. 

In summary, we proposed theoretically and confirmed in experiment that the quasiparti- 
cle transport in a superconducting SET transistor in the vicinity of the tunneling threshold 
can be described as resonant tunneling through a macroscopic charge state of the central 
electrode of the transistor. The maximal differential conductance associated with this pro- 
cess is e 2 /27rh, while the width of the resonance is determined by the lifetime broadening 
of the charge states of the transistor. For gate voltages away from the resonance we ob- 
served very clearly the co-tunneling current which exhibits linear (in contrast to cubic of the 
normal-metal case) dependence on the bias voltage. 

We gratefully acknowledge financial support of the Academy of Finland and US AFOSR, 
and Mikko Leivo for help with the experiments. 
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FIGURES 



Calculated bias-voltage dependence of the differential conductance of a symmetric su- 
perconducing SET transistor with junction resistance R = 20h/e 2 . The curves are plotted 
for several values of the gate voltage, i.e., the charge Q induced on the central electrode 
of the transistor, that correspond to several charging energy barriers E for tunneling: 
E /(hA/Re 2 ) = 0; 1; 3; 6; 10. The induced charge Q can not be close to 0. The inset 
shows the equivalent circuit of the SET transistor. 

Measured (a) I — V characteristics, and (b) bias-voltage dependence of the differential 
conductance of sample 1 for several gate voltages. The traces shown with thick lines in (b) 
correspond to the I — V curves presented in (a). For clarity, the features due to the current 
peaks associated with the Cooper-pair tunneling that are visible in (a) have been omitted 
in (b). For discussion see text. 

Measured I — V characteristics of sample 2 for several gate voltages. The inset shows the 
differential conductance in the vicinity of the gap edge V = 4A/e for gate voltages which 
correspond to the two largest tunneling thresholds. The conductance jump at V = 4A/e is 
due to the co-tunneling. 
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TABLES 



sample 


R(Ml) 


E c (meV) 


A (meV) 


a 




Go^theory 




1 


206 


0.35 


0.22 


0.08 


3.1 -10~ 3 


1.6 -10~ 3 


.5 


2 


152 


0.15 


0.21 


0.24 


0.032 


0.014 


0.9 


3 


65 


0.15 


0.20 


0.55 


0.096 


0.086 


1.6 


4 


52 


0.08 


0.23 


1.44 






4 .0 



TABLE I. Parameters of the four studied SET transistors. Conductances in the last three 



columns are shown in units of e 2 /2irh. 
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